Effective Operators for Double-Beta Decay 
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We use a solvable model to examine double-beta decay, focusing on the neutrinoless mode. After 
examining the ways in which the neutrino propagator affects the corresponding matrix element, we 
address the problem of finite model-space size in shell-model calculations by projecting our exact 
wave functions onto a smaller subspace. We then test both traditional and more recent prescriptions 
for constructing effective operators in small model spaces, concluding that the usual treatment of 
double-beta-decay operators in realistic calculations is unable to fully account for the neglected parts 
of the model space. We also test the quality of the Quasiparticle Random Phase Approximation and 
examine a recent proposal within that framework to use two-neutrino decay to fix parameters in 
the Ifamiltonian. The procedure eliminates the dependence of neutrinoless decay on some unfixed 
parameters and reduces the dependence on model-space size, though it doesn't eliminate the latter 
completely. 



I. INTRODUCTION 

Next-generation double-beta-decay experiments will 
provide data on the mass of the lightest neutrino 0] . But 
to interpret the data accurately we will need to better 
understand the nuclear physics that influences the decay 
rate. The problem of calculating double-beta-decay ma- 
trix elements has occupied nuclear theorists for a long 
time. The latest in a long history of neutrinoless-decay 
calculations are thought to approximate exact results to 
within factors of 3 or 4, but that estimate is not much 
more than an informed guess. In short, the calculations 
need to be improved. 

Almost all current calculations are performed ei- 
ther in the Quasiparticle Random Phase Approximation 
(QRPA) or the shell model. The two methods differ in 
the types of correlations they take into account. The 
QRPA can include the cumulative effects of small num- 
bers of particles or holes far above or below the Fermi 
surface, but the correlations are of a specific type. The 
shell model is restricted to configurations with low ex- 
citation, but correlates them in all possible ways. The 
QRPA calculations indicate that nearly empty orbits rel- 
atively far from the Fermi surface play a significant role 
in the nuclear ground states, implying that the shell- 
model is leaving something out. We are not referring here 
to the short-range correlations coming from the nuclcon 
hard core, the representation of which requires very high- 
lying excitations. QRPA correlations are longer range, 
connected to orbits just 10 or 20 MeV away from the 
Fermi surface; the shell model usually cannot include 
even these. Within the levels it does include, however, 
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the shell model indicates that QRPA correlations are too 
simple. We need either to add the effects of more com- 
plicated correlations to the QRPA or to include effects 
of additional single-particle levels in the shell model (or 
perhaps do both). 

The way to correct for the finite shell-model space is 
to construct effective Hamiltonians and transition opera- 
tors. Unfortunately, the equations that determine effec- 
tive operators are hard to solve and perturbative approx- 
imations often do not work. Phenomenology is the usual 
substitute. Thus, the effective Hamiltonians computed 
within perturbation theory usually include only two-body 
terms that are then renormalized to fit ground-state ener- 
gies and other observables. Effective transition operators 
are often taken to have the same form as the bare oper- 
ators but are multiplied by effective coupling constants. 
Such prescriptions can work well when there are lots of 
data to which parameters can be fit, but for the two-body 
operator governing neutrinoless double-beta decay, which 
has not yet been observed at all, the use of the bare opera- 
tor with or without simple multiplicative renormalization 
is suspect. Yet such is the state of the art in shell-model 
calculations of double-beta decay. To improve these cal- 
culations, we should focus on a better determination of 
the effective decay operator. 

Here we use a solvable model based on the algebra 
50(5) to test some of the assumptions that go into cal- 
culations of double-beta decay in ^^Ge, one of the most 
promising isotopes for experiment. We diagonalize the 
model Hamiltonian and calculate the rates of both two- 
neutrino (2j/) and neutrinoless (Oi') decay in a space con- 
sisting of a degenerate pfgg/2 valence shell, and a higher- 
energy degenerate scLg^ /2 shell. In this quite realistic 
model space (albeit with unrealistically simple correla- 
tions), we can compare the behavior of the two kinds of 
decays, which depend differently on radial wave func- 
tions, as we vary couplings in the Hamiltonian. We 



2 



can also project the full wave function onto states in 
which all nucleons are in the lower shell and test pro- 
cedures for constructing effective interactions and opera- 
tors. Though technically we are calculating Fermi decay, 
we allow isospin to be violated so that our calculation 
mimics many of the aspects of Gamow- Teller decay. We 
find, for example, that an effective two-body operator 
that reproduces decay rates exactly in nuclei with just a 
few valence nucleons, while better than the bare operator, 
is far from perfect when applied to a nucleus like ^^Ge 
that has many valence particles. The more phenomeno- 
logical prescription of multiplying the bare operator by 
an effective coupling requires different renormalization 
factors for 2v and Gv decay. Such results, while perhaps 
unsurprising, contradict common practice. 

We also examine some related ideas. A recently pro- 
posed technique for constructing effective operators non- 
perturbatively through a Lanczos approximation to the 
many-body Green's function usually proves efficient but 
appears to offer little advantage over the Lanczos meth- 
ods already used in shell-model calculations for double- 
beta decay. A recent proposal for adjusting the Hamil- 
tonian in the QRPA to fit measured 2v decay, though it 
eliminates the dependence on some of the parameters in 
the Hamiltonian and reduces the dependence on model- 
space size, does not work as well in this model as in more 
realistic calculations. The QRPA itself, ironically, repro- 
duces exact solutions extremely accurately here. 

The rest of this paper is divided as follows: In section 
^we discuss the model. Though it has been used before 
in 2v decay 0, ours is the first appHcation to Qv decay, 
the operator for which is not built from SO (5) generators 
and therefore takes some effort to manipulate. Section 
mil tests common prescriptions for constructing effective 
operators in a truncated model space against exact solu- 
tions in the full moded space. In section Hvl we examine 
the proposal of Ref. 3] to renormalize the Hamiltonian 
in the QRPA. Sectionals a conclusion. The appendix 
contains explicit formulae for matrix elements of double- 
beta decay operators in the model. 



II. SO(5) MODEL 

We consider the nuclei ^^Ge and ^^Se in a "small" shell- 
model space containing 12 protons and 24 neutrons (for 
Ge) in a degenerate pfgg/2 valence shell, and in a "large" 
space containing the small space together with all possi- 
ble excitations of the nucleons into a degenerate version 
of the next shell (sc?(77/2), which is an energy e above the 
pfgg/2 shell. We usually use one of two different values 
of e so that either about 20% of the particles are in the 
upper levels, or about 6% . Both values of the splitting, 
probably, are unrealistically small for the kind of long- 
range correlations we have in mind, but our correlations 
are unrealistically simple and we compensate in part by 
lowering e (for one test we will use a very large e that 
results in 1% of the particles in the upper levels). The 



single-particle wave functions come from a harmonic os- 
cillator with length parameter 6 = 2.12 fm. 

The Hamiltonian, constructed from the generators of 
SO{5) X SO{5) (one 5*0(5) for each degenerate set of 
levels), is 



H = eN2 (1) 

2 

— G ^ ^ {^pp^pp ~^ ^nn^nn dpp'^pn'^pn ^ 9phTa ' Tf,) , 

a,b=l 

where a, 6 = 1,2 label the shells (lower and upper), e is 
the energy difference between the shells, A'2 is the number 
operator for the upper shell, Ta is total isospin operator 
for shell a, and 

PP " 2 / , Jay^a^am 
'-'nn — 2 / .Jara'^aJo 

4n - ^Ei-t^-'^ilo- (2) 

Here tt^ (z^^) creates a proton (neutron) in level a with 
angular momentum ja, j = y/2j + 1, and the square 
brackets indicate angular-momentum coupling. The al- 
gebra of 50(5) contains the three pair-creation operators 
above (for a given set a of levels) , three corresponding de- 
struction operators, the three components of the isospin 
Ta, and the number operator Na- Since H contains only 
generators of 50(5) x 50(5), its lowest lying eigenstates 
will consist of configurations in which the nucleons are 
entirely bound in isovector 5 pairs of the type in Eq. (jj)). 

The Hamiltonian in Eq. is not the most complicated 
that could be built from generators of 50(5) x 50(5). 
One could, for instance, assign different strengths to the 
pairing forces in the two levels, or make the proton pair- 
ing stronger than the neutron pairing. Even without the 
last complication, the Hamiltonian violates isospin (T) 
conservation unless gpp = 1. In our application of this 
model, we will make an analogy between double Fermi 
decay, the 21^ version of which vanishes when isospin is 
exactly conserved, and double-Gamow- Teller decay in a 
more general model, which vanishes when Wigner SU (4) 
is exactly conserved. Of course, 51/ (4) is more badly 
violated than isospin, which is nearly a good quantum 
number, but in our simplified model, we will let isospin 
be badly violated to mock up the violation of 5'[/(4) in 
more realistic calculations'^. When T is violated, the 
space of fully-paired states in the two sets of levels pfgg/2 



SU{A) is violated primarily at the mean-field level, while in our 
calculations isospin is violated only in the two-body interaction. 
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and sdg'j/2 has dimension 1042 for ''^Ge and 1347 for the 
daughter nucleus ''^Se. 

We will calculate the ground-state to ground-state 
transition matrix elements of two operators: the 2u 
double-Fermi operator and the double-Fermi operator 
(both obtained from the closure approximation). These 
operators have the form 



(3) 



where i,j refer to particles and r"*" changes a neutron 
into a proton. The space/spin parts of the operators are 



and 
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(4) 



(5) 



The label k in Eq. refers to the kind of decay. The 
denominator in Eq. reflects the propagation of a vir- 
tual neutrino inside the nucleus. We have neglected the 
average nuclear excitation and short-range correlations 
in writing Oq^. We will also evaluate an operator that 
includes those effects: 



/(In-r,|) 



n - ri 



where 0,0 



(6) 



(7) 



with 7i — 1.1 fm^^, 72 — 0.68 fm^^, and E taken some- 
what arbitrarily to be 15 MeV. 

The fully-paired basis states can be labelled 
\Afi,Ti,Mi;M2,T2,M2), where Afi refers to the 
number of pairs in level-set 1 (the pfgg/2 shell), Ti 
to the isospin of those Afi pairs (with Afi — Ti even). 
Ml to the isospin projection of those pairs, etc., and 
Ml + M2 = 1/2(Z - N). The Hamiltonian, Eq. (HJ, and 
21^ decay operator, Eq. Q), are products of generators 
and one can evaluate their matrix elements, which 
depend only on the quantum numbers of the pairs, 
following, eg., Ref. But the Ov operators, Eqs. 

(O and lO, contain position-dependent factors that 
do not belong to the algebra. Their matrix elements 
will depend on the wave functions of the single-particle 
states that make up our basis and require more effort 
to evaluate. We have computed them by applying the 
generalized Wigner-Eckart theorem 0] , which for matrix 
elements of an operator M between fully paired states 
in SO{5) is: 

x{{na,0)\\M^'^"'^'^\ma,0)) {T^M^;ToMo\TaMa) 

x{{na,0)KTa-A^l,^2WoTo\ma,0)MaTa) (8) 



Here the extra "quantum number" fla, omitted in the 
labeling scheme discussed in the previous paragraph, is 
the half the total degeneracy of the level-set a, 
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(9) 



and labels the representations {ui = fla , 0^2 = 0) of 
50(5) in which all particles are fully paired. More gen- 
eral representations, such as that characterizing the op- 
erator in Eq. (|HJ, require two nonzero labels uji and uj2- 
The first factor on the right hand side of that equation is 
a reduced matrix element that depends only on fla and 
the operator quantum numbers lui and UJ2 ■ AH the depen- 
dence on the initial and final number of particles in the 
system, the initial and final isospin and isospin projec- 
tion, and the isospin and particle-number quantum num- 
bers of the operator appear in the double-barred 50(5) 
"reduced Clebsch-Gordan" coefficient (which is indepen- 
dent of the isospin-projection quantum numbers) and an 
ordinary isospin Clebsch-Gordan coefficient. 

The work of of Ref. 7] allows us to decompose 
the operators Mqi, and A^q^ into operators with good 
50(5) X 50(5) quantum numbers, lists some of the 
50(5) Clebsch-Gordan coefficients we need, and contains 
enough information to allow us to construct (laboriously) 
those it doesn't list: the coefficients for isospin-changing 
components of operators belonging to the 50(5) repre- 
sentations (2, 0) and (2, 2). Here we show expressions for 
matrix elements of double-charge-changing two-body op- 
erators in the "small space" — all particles in the degen- 
erate pfgg /2 levels — the states of which we write in the 
form I A/", T, M) , where all quantum numbers refer to the 
first set of levels (a = 1) because the second set is empty. 
We focus on these matrix elements because they are the 
most complicated and because they will determine the 
behavior of the effective decay operators when we trun- 
cate the wave function to the small space (other matrix 
elements needed for the full calculation are in the ap- 
pendix) . Defining summed particle-particle and particle- 
hole-like matrix elements of an arbitrary two-body oper- 
ator A^K that changes charge by two units (acting, for 
generality of notation, in either of the two sets of levels) 
as 



where, as before, a and b refer to a degenerate sets of 
levels and the matrix elements are antisymmetrized, we 
have when all the particles are in the lowest set 
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x(T±2 M + 2;l - l|r> - 1 M + 1){T M; 1 1|T> - 1 Af + 1) (2J^p''[1,1] - (2^2 - 1]) , 

and 



(A/", T, Af + 2|A/{«,|A/', T, M) = — ^^^^^"^^ (T Af; 2 2|r Af + 2) 

^ I -I ' ' / 4^6(2T - l)(2r + 3)r!(r! - l){2n + 1) ^ ' ' 



:{[AA(AA-2J7)-3(2f7 + l) + (4f7 + l)T(r+l)] 2J^p''[1,1] (12) 

- [{2VL - i)N{N - 2Q) - 3n{2n + 1) + (20 + 3)r((T + 1)] jc-pp[i^ 1] I . 



Here f2 = fii, and T> (r<) is the largest (smallest) of the 
isospins in the bra and ket of Eq. p2(l . The eigenstates 
of the Hamiltonian Eq. ^ have good particle number Af 
and isospin projection Af , but are mixtures of states with 
different isospins T. The same is true of the projection 
of the wave functions onto the small space, so both the 
equations above are important. 

The most important fact about Eqs. pi|) and p2|) is 
that they depend on the two-body matrix elements of 
the operator through only two quantities, ^^^^[1,1] and 
JFP''[1,1], the particle-particle and particle-hole (with 
some factors of j removed) matrix elements summed over 
all orbits in the lower shell. The action of any two-body 
operator that changes charge by two units on fully paired 
states in that shell can therefore be constructed by spec- 
ifying just two parameters. This simplicity reflects the 
nature of the underlying correlations. The interaction 
cannot create particle-hole excitations involving different 
orbits. 



III. MATRIX ELEMENTS, TRUNCATION OF 
WAVE FUNCTIONS, AND EFFECTIVE 
OPERATORS 

Before addressing the issues of effective decay opera- 
tors in the small space, we examine the behavior of the 
matrix elements 



M^ = {Of\M^\0^ 



(13) 



of the operators Ai2iy, AAqu, and AA'q^j connecting the 
ground states of ^^Ge and ^^Se. The matrix elements 
are plotted in Fig. ^ for fixed G, e = lOG (resulting 
in about 20% of the nucleons occupying the upper set of 
levels), and two values of gph. We vary gpp, which mea- 
sures the strength of proton-neutron pairing (the graphs 
look qualitatively the same when we change e to 20G and 
about 6% of the particles are in the upper set.) To com- 
pare the behavior of the three matrix elements with gpp, 
we have scaled the matrix elements so that they are 
equal to the 21/ elements at gpp = 0. At gpp = 1 isospin 



is conserved and the 2v matrix element crosses zero. In 
realistic calculations the analog of gpp in the J'^ = 1+ 
channel is close to 1; indeed, early QRPA calculations 
gave matrix elements that were much too large precisely 
because they neglected the particle-particle interaction. 
The matrix elements cross zero only for larger gpp, 
much larger when gph = 1, and are therefore less sensi- 
tive than the 21/ elements to that parameter when it is 
at a realistic value. This reduced sensitivity is another 
aspect of our model that mirrors more complicated cal- 
culations. The two versions of the Qv matrix element are 
nearly proportional; the short-range correlations and av- 
erage excitation energy shrink the matrix element by a 
factor of about 1.5, almost independently of gpp. 



in 



o 

GO 



10 


-10 



































— M'q^ (x 3.6 fm'') 

- - (X 2.3 fm"') 








1 


























— (X 2.5 fm"') 








- - (X 1.7 fm ') 


gph=l 






1 







0.5 



FIG. 1: Double-beta-decay matrix elements versus the 
strength of neutron-proton pairing g^^ for two values of the 
particle-hole coupling g^y^. The matrix elements of the two 
operators are scaled to that of the 2v operator at g^^ = 0. 
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The different behavior of the 2v and Qv matrix ele- 
ments with gpp is explored in Fig. Q, which shows the 
pair-separation transition density 



P(r)^(0/|^<5(r-|n-r,|)r+r+|0.) , 



(14) 



(r is the magnitude of the internucleon separation) at 
several values of gpp. The density P{r) is defined so that 



P{r)dr ^M2u 
P{r) 



-dr = M. 



Of 



(15) 
(16) 



etc. The figure shows that as gpp increases P(r) decreases 
and then changes sign at large r, while remaining positive 
at small r. At gpp = 1 the integral of P{r) is zero but the 
1 /r- weighted integral is still positive because small values 
of r are most important. Thus, the matrix elements of 
Mou and M'f)^ cross zero at larger gpp than that of Al2i/- 
The reason for that P{r) is most affected by gpp at 
large r is the following: At gpp — only one pair can par- 
ticipate in the ground-state to ground-state decay; since 
there are no correlated proton-neutron pairs, the decay 
operator turns two neutrons within a given pair into two 
protons. As gpp becomes larger the ground states contain 
neutron-proton pairs, and double-beta decay can pro- 
ceed, e.g, through the transformation of one neutron in 
a neutron-neutron pair and one in a (separate) proton- 
neutron pair. These two neutrons are farther apart on 
average than two in the same correlated pair, so the tran- 
sition density changes more at large r than at small r, 
where it stays positive. 
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FIG. 2: The transition density P{r) (see eqn. dJ) as a 
function of internucleon separation r for several values of gpp, 
with e = lOG and gph = 0. 

We now turn to the question of effective operators in 
the small space (all nucleons in the pfgg/2 shell). The 



Bloch-Horowitz equation 



H,s{E)^PHP + PHQ 



1 



E-QH 



QHP 



(17) 



and approximations thereto are designed to construct an 
effective Hamiltonian _ffcff that obeys 



(18) 



where H is the full Hamiltonian, j^Pfe) is an eigenstate 
with energy E^., P projects onto the subspace in which 
all particles are in the pfgg/2 shell, and Q = 1 — P. Once 
the projections of the eigenstates are known, effective 
transition operators can be constructed to satisfy 



where 



PI*) 



(19) 



(20) 



is the normalized component of \^) that has all the par- 
ticles in the lower set of levels. We don't actually have to 
know HcS to construct Al^*^; we just need the full wave 
functions \^k) and their normalized projections \^k)- 

Effective operators must in general be sums of 1-body, 
2-body, ... up to N-body terms. It is common practice, 
however, to attempt to incorporate the effects of the ex- 
cluded space in an effective operator of the same rank 
(number of bodies) as the original operator. Thus effec- 
tive electromagnetic transition operators are often taken 
to be one-body but of a more general form than the orig- 
inal operator, with "level-dependent effective charges". 
In calculations of 2v double-beta decay, people often 
just multiply the bare operator by a single effective cou- 
pling, usually 1/(1.25^), to take into account the observed 
quenching of spin-dependent transitions. When calculat- 
ing Qv decay they either use the bare operator or the 
same quenched coupling as for 2v decay. Nothing beyond 
a simple multiplicative renormalization has ever been at- 
tempted. 

Although by definition there always exists a multiplica- 
tive factor that will work, the trick is to know what it is. 
Figure 13 shows that the appropriate factor is different in 
Ov decay than in 2v decay. At gpp — 0.9, a value yielding 
a reasonable approximation to the real world within our 
model, the renormalization factor for Qv decay must be 
1.8 when e = lOG and 1.5 when e = 20G {gph = 0). For 
2v decay, these numbers are 1.3 and 1.1 (adding the par- 
ticle hole force changes the numbers only a little) . These 
results measure the quality both of calculations that use 
bare operators and of those that attempt to incorporate 
the effects of neglected single-particle levels by renormal- 
izing the 2v and operators by the same factor. 

Figure 21 which shows P{r) at gpp = 0.5 for the fuU 
ground states |0) and their normalized small-space pro- 
jections |0), sheds some light on the different renormaliza- 
tion factors. Interestingly, the transition density spreads 
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FIG. 3: 2v (top) and Qu (bottom) matrix elements versus the 
strength of neutron-proton pairing gpp for Qph = 0. The solid 
lines are the results in the full models space and the dashed 
lines are the results obtained by sandwiching the bare opera- 
tor between normalized projections of the wave functions onto 
the small model space. 



out when projected onto the small space, so that most of 
the change is at small r and, as a result, the renormaliza- 
tion factor for the Qv transition is larger than for the 2v 
transition. One can understand the spreading by looking 
at the structure of the correlated pairs in Eq. ||2Jl. The 
one-body piece of the mean-square internucleon distance 
((ri — r^) ■ (ri — r2)) is uniformly reduced by the cross 
term — 2(ri • r2). In the pfgg/2 shell r connects only the 
fi/2 and ggf2 levels. Adding the upper shell introduces 
many more contributions to (ri • r2) and so reduces the 
mean-square radius below its value in one shell. Alterna- 
tively (and more simply), the Oz/ matrix element is more 
sensitive to occupation of the upper set of levels because 
the neutrino propagator in Eq. jSJ allows the Ov operator 
to connect those levels with the more highly occupied lev- 
els in the lower shell. The operator, which lacks radial 
dependence, does not connect levels from different shells. 

If simple renormalization by a constant factor is dan- 
gerous, what about applying ideas used in electromag- 
netic physics to obtain a more general two-body opera- 
tor? The usual practice there is to assume N nucleons 
in the valence shell when constructing an effective iV- 
body operator, and then to apply that operator when 
there are larger numbers of valence particles. We can 
implement that procedure easily and exactly here, with 
N = 2. The ground-state to ground-state transition for 
two valence neutrons in the small space decaying to two 
protons doesn't depend on the parameter ^^'l.g of the ef- 
fective operator in Eq. 119|) (i.e. the coefficient of .^-"^^[1, 1] 
in Eq. (|12|l vanishes when JV — 1 and T — 1 because the 
particle-hole piece of the operator doesn't act on nucleons 
within the same pair), and so determines the parameter 
„. The excited states in both the initial and final 

At, en 

nuclei have seniority two and are unaffected by the in- 
teraction; their structure plus that of the paired ground 



FIG. 4: The transition density P(r) as a function of internu- 
cleon separation, at gpp = 0.5, for the full ground-state wave 
functions (solid line) and their normalized projections onto 
the small space (dashed line). 



state determine ^'l^g in a simple way. More precisely, 
the effective particle-particle and particle-hole parame- 
ters for the operator M'^ are given by 



{Of\M^\0^ 



(0/|Xk|0,; 



and 



2n 



;,cff 



(21) 



(22) 



where the wave functions \0i) (|0/)) arc the full two- 
neutron (two-proton) ground-state wave functions, the 
barred wave functions are their normalized projections 
onto the small space, and 17 = ili characterizes the 
pfgg/2 shell. With these identifications (the JF[l,l]'s 
are "bare" coupling constants), the operator M'/^ repro- 
duces all transitions involving the lowest 17(217—1) states 
in each nucleus. It is what one would obtain by carry- 
ing the usual diagrammatic perturbation theory for the 
two-valence-nucleon system to all orders. 

Of course we really ought to be determining the ef- 
fective operator for ^^Gc — >^^Se, rather than the two- 
valence-nucleon transition "'■^Ca — >^^Ti. That, however, 
would amount to solving the problem exactly. But one 
might decide that it would be better (and feasible) to fo- 
cus on ground-state to ground-state transitions, in the 
two-nucleon system and, say, in the four-nucleon sys- 
Ti. Reproducing these two transitions 



tern 44 Ca 



would determine the two parameters 3-^^^^ and ^^^^g, 
and would have the advantage of incorporating neutron- 
proton pairing, which we know plays an important role 
in the heavier systems when gpp is turned on, into the 
effective operator (through the two-proton two-neutron 
state in '*''Ti). 

We display the results of both these prescriptions for 
the operator Mq^ in Fig. El Both procedures improve 
on the bare operator, particularly when e = 20G, where 
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FIG. 5: Matrix elements of Ov effective operators as a func- 
tion of Qpp. The dashed line corresponds to the bare oper- 
ator in the small model space, the dotted line to the effec- 
tive operator that reproduces transitions in the two-valence- 
nucleon system, and the dot-dashed line to the operator that 
reproduces ground-state transitions in two- and four-valence- 
nucleon systems. The solid line is the exact result in the full 
model space. 



they make up about half the difference between the bare 
and full results, but neither is perfect^. As expected, 
the fit to transitions involving two and four nucleons 
works better than the operator constructed entirely in 
the two-valence-nucleon system, but not remarkably so. 
We tried determining the effective operator entirely from 
four-valence-nucleon systems, but the improvement was 
minimal. The difference between the exact and approx- 
imate results is a measure of the amount higher-body 
effective operators contribute to the transitions. 

One reason that two-valence-nucleon systems are em- 
phasized in the traditional treatment of effective opera- 
tors is that diagrammatic perturbation theory becomes 
complicated when more particles are considered. Refer- 
ence proposed a nonperturbative algorithm, based on 
the Lanczos representation of Green's functions as con- 
tinued fractions [lol | , that can be used for any number of 
particles. For every state \Af, T, M) in the small space 
and some guess for the ground state energy E (in both 
the initial and final nucleus), we construct the approxi- 
mation to the large space vector {E — QH)^^\Af,T,M) 
and Hcs{E), and then diagonalize the latter in the small 
space. We then redefine E to be the lowest eigen-energy 
and repeat the process until E doesn't change. Finally we 
construct the Lanczos approximation to the full ground 



^ For 2v decay the bare operator gives a better approximation than 
for Ov decay (see Fig. UJ, but the effective-operator procedures 
improve on it only slightly. 
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where is a computable normalization factor, and cal- 
culate the transition matrix elements. The results for 
different numbers of Lanczos vectors appear in Fig. |S1 
where e — 20G. Only 8 or 9 such vectors are required 
to get excellent results, but the iteration method fails to 
converge if e is too small (e.g. lOG.). Furthermore, it 
is hard in our context to see the advantage of this pro- 
cedure over simply doing a Lanczos diagonalization to 
obtain the ground-state eigenvectors in the large space 
with the "starting vector" determined, for example, by 
diagonalizing H in the small space. Such techniques are 
already exploited in shell- model calculations, and it's not 
clear that replacing them with the Bloch-Horowitz equa- 
tion would be useful. Reference pj| proposed a new per- 
turbative approach to solving that equation, but it has 
been worked out only in very light systems. 
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FIG. 6: Approximations to the Oi/ matrix element, for e = 
20G, from a Lanczos representation of the Bloch Horowitz 
equation, as described in Ref. Results are shown when the 
representation is truncated at different numbers of Lanczos 
vectors, and compared to the exact results. 



IV. QRPA 

The QRPA is the most commonly used method for cal- 
culating double-beta decay rates. How accurate is it? If 
it is not so accurate, can its predictions for Ov decay be 
improved by fixing certain parameters in the Hamilto- 
nian to reproduce 2v decay, as suggested in Ref 0? We 
examine both questions here. 

We carry out the QRPA in the usual way (see, e.g., Ref 
5]). We first solve the BCS equations in the initial and 
final nuclei, obtaining the amplitudes u^^ , , , 
and , . The corresponding occupation 



probabilities agree very well with the exact occupation 
probabilities for all values of the energy splitting e. Next, 
we solve the usual QRPA equations of motion for both 
nuclei and use the expressions in Ref. Q to evaluate the 
transition matrix elements. (The only difference is that 
in that paper the matrix element was obtained by aver- 
aging results from a calculation purely within the initial 
nucleus and from one within the final nucleus. Here we 
use both sets of u's and f 's in the same expression, in- 
serting an overlap between intermediate states generated 
by the QRPA from the initial and final ground states.) 
Since we evaluate only the closure matrix elements we can 
compare with the exact solutions. The dependence on 
the level splitting e resides in the BCS amplitudes u, v as 
well as in the QRPA amplitudes X, Y; the latter also de- 
pend on the interaction strength gpp. The QRPA matrix 
elements are always evaluated by expanding the transi- 
tion operator in proton-particle/neutron- hole multipole- 
mutipole form. Since for 2iy decay the transition operator 
is just 2tiT2, only the 0+ contribution (in real nuclei, 
the 1"*" contribution) exists. The neutrino propagator 
causes all multipoles contribute to Ov decay. 

In our model the interaction acts only in 0^ intermedi- 
ate states. For all other multipoles Ov expressions reduce 
to their BCS form. The final result can be expressed as 

(0/|A<,|0,)qrpa = (0+ part) (24) 

~. c J'^SaR 
X JaJl30jfi ^- , 

where the matrix elements with rounded brackets are un- 
symmetrized (as usual in the RPA) and "0+ part" refers 
to the matrix element when only the 0^ multipole in the 
particle-hole decomposition of the operator is included. 
For 21^ decay the rest of the expresson vanishes, but for 
Ov decay there is a simple contribution from the other 
multipoles that is independent of the proton-neutron in- 
teraction. 

The exact and QRPA results for both double-beta de- 
cay modes are compared in Fig. [7| The agreement is 
very good over the whole interval gpp = 0—1. Note 
that QRPA equations have no solution (the solutions are 
said to "collapse" ) for gpp only slightly above unity. The 
agreement is equally good for all values of the splitting 
e. Previous comparisons between QRPA and exact re- 
sults for M21, in this model P| have not found such good 
agreement. The better agreement here is related to the 
larger degeneracies and the larger number of particles 
we use to represent shell- model calculations in ^^Ge. 

Reference suggested a procedure for reducing the 
dependence of the Ov rate on the number of single- 
particle states used in the QRPA. Simply put, the 
method is this: assuming that the strength of the 
particle-particle interaction, gpp, is the most important 
parameter, one evaluates both matrix elements Mqi, and 
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FIG. 7: Comparison of the exact and QRPA values for the 
2v matrix element (top) and Ov matrix element (bottom) for 
the indicated values of the spacing e between the degenerate 
fpgg/2 and sdgY/2 shells. 



as a function of that parameter. One can then in- 
vert the relation between and gpp and express Mq^, 
directly in terms of M2u- Reference 3] did this within 
the QRPA and in a modified version, the Renormalizcd 
QRPA [13, for 76Ge, i°°Mo, "Oxe, and ^^exe. It turned 
out that the effects of single-particle levels far away from 
the Fermi surface, while noticeable for Mq^, and M21, con- 
sidered separately, became negligible when Mqi/ was plot- 
ted against M2i^- In other words, knowing M21, (from 
experiment) one could obtain the unknown Mq^, and its 
value was the same whether the far away single-particle 
states were included or not in the QRPA, as long as at 
least two oscillator shells (more than in our small space) 
were included. The question remains whether the QRPA 
result, even if insensitive to the model space, is correct. 

When we carry out an analogous procedure in the 
50(5) X 50(5) model, we find that Moi, and M21, are 
related as shown in Fig. |S| While the relation between 
them is a simple one, it is not independent of e, as it 
would be if the ideas of Ref. 0| played out perfectly. In- 
stead, the offset in Mq^ depends on e (but is independent 
of gpp). With increasing e the lines in the figure, for ob- 
vious reasons, come closer together. While fixing gpp to 
the exact (large-space) value of M2,y improves the small- 
space prediction of Afoiy, it always leaves it smaller than 
the exact result. 

In the realistic calculations the strength of the pair- 
ing interaction was adjusted so that the pairing gap was 
the same in the "large" and "small" spaces. We do the 
same a.t e — 20G; the corresponding line moves closer to 
the e = 00 limit, but is still significantly different. Fi- 
nally, in Ref. Q the proton-proton and neutron-neutron 
pairing constants were adjusted separately. With a typ- 
ical value G ^ 22/A MeV and huj = 41 /A^/^ MeV, the 
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value of e corresponding to Ref. |3| is at least 40G. With 
that value and a charge-dependent pairing renormaliza- 
tion Gnn/Gpp ^ 1-2 here, the line moves even closer to 
its limiting value. 

The matrix elements in Fig. [S] were obtained without 
any particlc-holc force, i.e. for gph — 0. When the cal- 
culation is repeated with gph = 1 both Mq,^ and M21, 
are significantly reduced (see Fig. ^| . When the param- 
eter gpp is eliminated, however, and Mq^ is expressed 
as a function of M2j/, the resulting lines are essentially 
identical to those shown in Fig. |H1 In other words, the 
procedure of Ref. Q eliminates the dependence of Mqi, 
on the parameter gph, at least in this model. This impor- 
tant conclusion should be tested in more realistic QRPA 
calculations. 
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FIG. 8: Dependence of the matrix element Mou{gpp) on the 
matrix element M2v{gpp) for the indicated values of the level 
spacing e. See text for further explanation. 



truncation, 2) familiar procedures for finding effective 
two-body decay operators to correct for truncation help 
matters but leave substantial room for improvement, and 
3) The procedure of Ref. 3] for eliminating model-space 
dependence in the QRPA helps but doesn't work as well 
in our model as in realistic calculations. 

Several issues still need to be investigated. How ac- 
curate are the QRPA results of Ref. 3] once the renor- 
malization procedure is applied? Ironically, although the 
QRPA itself is very accurate in our model, almost cer- 
tainly more accurate than in reality because the model 
has very collective modes (the pairs), the renormaliza- 
tion procedure leaves more residual model-space depen- 
dence here than in realistic calculations. We therefore 
can't really address the question in our model. Another 
issue we haven't addressed is the effects of very-high mo- 
mentum states on double-beta-decay operators. Such ef- 
fects are usually simulated by the short-range-correlation 
function / in Eq. Q, but it's not clear how accurate a 
result is achieved that way. We are currently examining 
two-body short-range correlations by summing ladder di- 
agrams similar to those that make up the Brueckner G- 
matrix. 

Properly treating more-body effects will be harder, 
both for short-range correlations and the longer-range 
correlations we investigated here. Finding a good ap- 
proximation scheme to incorporate them into shell model 
calculations is, in our view, the most important problem 
in the theoretical treatment of double-beta decay. Be- 
cause only the ground states need to be calculated, much 
of conventional effective-operator theory is superfluous 
and the biggest payoff may come from pushing Lanczos- 
diagonalization methods to the largest model spaces pos- 
sible, leaving only high-momentum correlations to be ab- 
sorbed into operators. 



Why does the prescription of Ref. 3] fail to fully elim- 
inate model-space dependence in this model? There are 
several possibilities. As we pointed out earlier, the mag- 
nitude of the level splitting e used here is unrealistically 
small. Using larger e moves the large-space lines closer 
to those of the small space. Also, the interaction here af- 
fects only the 0+ multipole, while in realistic calculations 
all multipoles are affected to some degree. In any event, 
even with its shortcomings the procedure is interesting 
enough here that we find Afoi/ = a x + h indepen- 
dently of the values of both gpp and gph, with h decreasing 
as the level splitting increases. Combined with the ac- 
curacy of the QRPA in our model, this implies that the 
procedure gives a lower limit for the exact value of Mq^, 
provided that the Hamiltonian is corrected to reproduce 
the exact value of M2u- 

V. CONCLUSIONS 

Our most important findings are that 1) calculated 2v 
and Qv decay rates are affected differently by model-space 



APPENDIX A 

Here we supplement Eqs. (jllf) and (|12|l with rela- 
tions necessary to evaluate matrix elements of charge- 
changing two-body operators between states of the form 
\Ni,Ti,Mi;N2,T2,M2), corresponding to particles in 
both sets of levels. Eqs. Hll|l and 112|l can be used with 
M ^Mi {M2), T = Ti (T2), M = Ml (M2) and the [1,1] 
([2,2]) rs whenever A/'a (Mi), T2 (Ti), and M2 (Mi) 
are the same in the bra as in the ket. When those 
isospin quantum numbers change, we also need [with 
M[ = Ml + 1 and M2 = M2 + 1 everywhere below): 

{Ni,T[, M[;Af2, T!2,M^\M,Wi,Ti, Mi; T2, M2) = 

1. . Vri(ri + i)r2(T2 + i) 

-2^Tun5T.,T- (Al) 

x(Ti Mi;l llTi M[){T2 M2;l l|r2 M^) .ff [12] , 

where the delta function arises from the condition that 
M —T he even. When the operator moves pairs from one 
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set of levels to another, finally, we need and 



(M + 1, Ti, M[-M2- 1, % M^l A4«|M , Ti, Mi; M, T2, M2) 



1 



x(Ti Mi; 1 1|T^ M{)(T^ M^; 1 - 1|T2 M2) ^ril2] , V 2 1 

and 

- l,Ti',M{;AA2 + l,r^,M^|7W«|M,Ti,Mi;AA2,T2,M2) 



■^^^,-l.Ti.T^■^^l,.T■,.TL i^^) 



x(Ti' M{; 1 - l|Ti Mi)(r2 M2; 1 1|T^ M^) ^^[12] 
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